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Abstract 

We review recent progress in our understanding of clastic and in- 
clastic pion reactions on two nucleon systems from the point of view of 
effective field theory. The discussion includes 7rd scattering, jd — > 7r + nn, 
and NN — > NNir. At the end some remarks are made on strangeness 
production reactions like 7c? — > KAN and NN — > KAN. 

1 Introduction 

Even after several decades of research, the interactions and dynamics of strongly 
interacting few-nucleon systems are still not fully understood. Although phe- 
nomenological approaches are quite often very successful in describing certain 
sets of data, a coherent overall picture with clear connection to the fundamen- 
tal theory, QCD, is still lacking. 

The only way to a systematic and well controlled understanding of hadron 
physics is through the use of an effective field theory. The effective field the- 
ory (EFT) for the Standard Model at low energies is chiral perturbation the- 
ory (ChPT). This EFT already provided deep insights into strong interaction 
physics at low energies from systematic studys of the tttt [1] and the ttN [2, 3] 
system — reviewed in Ref. [4] — as well as the NN system [5-7]. Here we 
focus on the field of elastic and inelastic reactions on few nucleon system. 

A first step towards a systematic study of elastic and inelastic reactions on 
nuclei was taken by Weinberg already in 1992 [8]. He suggested that all that 
needs to be done is to convolute transition operators, calculated perturbatively 
in standard ChPT, with proper nuclear wave functions to account for the non- 
perturbative character of the few-nucleon systems. This procedure looks very 
similar to the so-called distorted wave born approximation used routinely in 
phenomenological calculations, but, in contrast to this, opens up the possibility 
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to use a power counting scheme. Within ChPT this idea was already applied to 
a large number of reactions like nd — ► ird [9], "fd — ► 7r°<i [10, 11], 7r 3 He^ 7r 3 He 
[12], 7r~(i — > 7nn [13], and 7c? — > ir + nn [14], where only the most recent 
references were given. We start our presentation with a brief description of 7rd 
scattering in sec. 2. 

Using standard ChPT especially means to use an expansion in inverse 
powers of Mjy — the nucleon mass. However, some pion-few-nucleon dia- 
grams employ few-body singularities that lead to contributions non-analytic 
in ttl^/Mn, with for the pion mass [15]. This is discussed in detail in 
sec. 3. There we show that the appearance of these contributions is linked to 
the Pauli principle operative while there is a pion in flight. In this context we 
discuss also the reaction jd — ► ir + nn, for in the mentioned sense this reaction 
is complementary to ird scattering. We show that within ChPT the existing 
data can be described to high accuracy and therefore the reaction qualifies as 
a tool to measure the nn scattering length. 

A problem was observed when the original scheme by Weinberg was applied 
to the reactions NN —> NNit [16,17]: the inclusion of potentially higher 
order corrections were large and lead to even larger disagreement between 
theory and experiment than found in earlier phenomenological studies [18]. 
For the reaction pp — > ppir® one loop diagrams that in the Weinberg counting 
appear only at NNLO where evaluated [19, 20] and they turned out to give 
even larger corrections putting into question the convergence of the whole 
series. However, already quite early the authors of Refs. [21, 22] stressed that 
an additional new scale enters for NN — > NNir that needs to be accounted 
for in the power counting. Since the two nucleons in the initial state need to 
have sufficiently high kinetic energy to put the pion in the final state on-shell, 
the initial momentum needs to be larger than 

Pthr = y/M N m n . (1) 

The proper way to include this scale was presented in Ref. [23] - for a 
recent review see Ref. [24]. As a result, pion p- waves are given by tree level 
diagrams up to NLO and the corresponding calculations showed satisfying 
agreement with the data. However, for pion s-waves loops appear already at 
NLO [23, 25]. In sec. 5 we discuss their effect on the reaction NN — > dir near 
threshold. In some detail we will compare the effective field theory result to 
that on phenomenological calculations. 

The central concept to be used in the construction of the transition oper- 
ators is that of reducibility, for it allows one to disentangle effects of the wave 
functions and those from the transition operators. As long as the operators 
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are energy independent, the scheme can be applied straight forwardly [26], 
however, as we will see below, for energy dependent interactions more care is 
necessary. This will also be subject of sec. 5. 

Once the reaction NN — > dir is understood within effective field theory 
one is in the position to also calculate the so-called dispersive and absorptive 
corrections to the ird scattering length. This calculation will be presented in 
section 7. 

When switching to systems with strangeness one immediately observes that 
the initial momentum needs to be quite large in any production reaction where 
there are two nucleons and no strangeness in the initial state. Strangeness 
conservation demands that at least two particles with strangeness in the final 
state. Therefore the mass produced is at least 

A = M A - M N + m K -m , 

where rrii denotes the mass present in the initial state in addition to the two 
nucleons, e.g., rrii ~ for ird — > KAN and rrii = for NN — ► NAK. There- 
fore, in the latter reaction the expansion parameter of the ChPT, namely the 
initial momentum in units of the nucleon mass, is 0.85, using the analog of 
Eq. (1). Also in the former reaction we are faced with an initial momentum 
0.6 in units of the nucleon mass — again not useful as an expansion parame- 
ter. Clearly, in those cases the chiral expansion as proposed above is no longer 
applicable. In sec. 9 we discuss, how for those large momentum transfer reac- 
tions the scattering lengths of the outgoing baryons can still be extracted in 
a controlled way using dispersion theory. 

We close with a brief summary and outlook. 

2 Remarks on the ird system 

We start our discussion with some remarks on the ird system. Pion-deuteron 
(ird) scattering near threshold plays an exceptional role in the quest for the 
isoscalar ttN scattering length a+, since the deuteron is an isoscalar target. 
Therefore one may write Re(a 7r ^) = 2a + + (few-body corrections) . The first 
term ~ a+ is simply generated from the impulse approximation (scattering 
off the proton and off the neutron; diagram (a) of Fig. 1) and is independent 
of the deuteron structure. Thus, if one is able to calculate the few-body 
corrections in a controlled way, ird scattering is a prime reaction to extract a + 
(most effectively in combination with an analysis of the high accuracy data on 
pionic hydrogen). In addition, already at threshold the ird scattering length 
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a) b) c) 

Figure 1: Typical diagrams that contribute to ird scattering. Shown is the 
one-body term (diagram (a)), the leading two-body correction (diagram (b)) 
and a possible short-ranged operator (diagram (c)). 

is a complex-valued quantity. It is therefore also important to gain a precise 
understanding of its imaginary part - this issue will be discussed in sec. 7. 
Recently the ird scattering length was measured to be [27] 

a^ p = (-26.1 ± 0.5 + i(6.3 ± 0.7)) x 10~ 3 m' 1 , 

where m n denotes the mass of the charged pion. In the near future a new 
measurement with a projected total uncertainty of 0.5% for the real part and 
4% for the imaginary part of the scattering length will be performed at PSI 
[28]. Clearly, performing calculations up to this accuracy poses a challenge to 
theory that several groups recently took up [29-34]. In addition, an interesting 
isospin violating effect in pionic deuterium was found, see [35]. For a review 
on older work we refer to Ref. [36]. 

A typical few-body correction to the ird scattering length is shown in di- 
agram (6) of Fig. 1. As we will see below, the contribution of this diagram 
largely exhausts the value of the ird scattering length not leaving much room 
for a contribution from a+, or stated differently, pointing at a small value of 
a+. Based on calculations within pion less EFT, it was claimed recently that 
this diagram is sensitive to the short range part of the deuteron wave func- 
tion [37,38]. As a consequence field theoretic consistency requires that at the 
same order there is to be a local operator to absorb this model dependence — 
the corresponding diagram is shown as diagram (c). Since this diagram comes 
with an a priori unknown strength not fixed by symmetries, ird scattering 
would be useless for the extraction of a+. However, systematic investiga- 
tions showed that, as soon as the pion exchange is included explicitly in the 
NN potential, diagram (6) can be evaluated in a controlled way [30,32-34]. 
Given this we assume from now on that short-ranged operators contributing 
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Figure 2: Formally leading few-body corrections to the ird scattering length. 

to the ttNN — > ttNN transition potential scale naturally. In other words, 
contribute with strength parameters of the order of one. Based on this we 
may estimate the contribution of diagram (c) of Fig. 1 relative to diagram (6) . 
This kind of analysis gives a relative suppression of the order 0(x 2 ), where 
X = m n /MN is the standard expansion parameter of ChPT with Mn (m n ) for 
the nucleon (pion) mass. This, together with the knowledge that diagram (a) 
largely exhausts the value of the ird scattering length, allows us to estimate 
the theoretical limit for the extraction of a+ from a measurement of the ird 
scattering length. We find 

Aa thco ~ 5 x 1(T 4 m- 1 . (2) 

To meet this theoretical limit we need to include in the calculation all contri- 
butions to the 7tg? scattering length lower than 0(x 2 )- 

Already in his original work, Weinberg discussed ird scattering at threshold 
as an illustrative example [8]. As usual, the leading contributions to the tran- 
sition operators are all those tree-level diagrams that can be constructed from 
the leading irN and irir Lagrangians. Those are shown in Fig. 2. Note that it 
is very important that the complete set of diagrams is considered, since, for 
example, both diagram (b) as well as diagram (c) are depending on the partic- 
ular choice made for the pion field. However, the sum of both is independent 
of this choice, as was first pointed out in Ref. [39]. 

Since all diagrams of Fig. 2 contribute to the same chiral order, naively 
one expects them to give similar contributions. However, explicit numerical 
evaluations showed, that diagram (a) exceeds the sum of the other two by 
about two orders of magnitude [8]. Can we understand this? One possible 
explanation could be that the sum of (6) and (c) is small because of significant 
cancellations between the two, probably due to the mechanism indicated at 
the end of the previous paragraph. Another possible explanation was given in 
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Figure 3: Typical pion loop contributions to trd scattering. Since the exchange 
of the two nucleons in the intermediate state (at the perpendicular lines) 
transforms one diagram into the other, the imaginary parts (and their analytic 
continuation) are linked by the Pauli principle, as described in the text. 

Ref. [9]: as a consequence of the small binding energy of the deuteron, e, the 
typical nucleon momentum inside the deuteron, 7, is also small, 7 ~ \/M^e, 
which turns out to be numerically about 1/3 of the pion mass. Since diagram 
(a) is proportional to the expectation value of l/q 2 , where q denotes the 
momentum transfer through the pion, and the sum of diagram (6) and (c) is 
proportional to the expectation value of q 2 /(q 2 + w^) 2 , we expect the ratio 
of the two contributions to be of the order of ( , -//m 7T ) 4: ~ 10 -2 , where \q \ was 
identified with the value of 7 defined above. Thus, the small binding energy of 
the deuteron seems to provide a natural explanation of the relative suppression 
of diagram (6) + (c) to (a). However, more systematic studies are necessary. 

Another important issue is the role of nucleon recoil contributions. All 
calculations mentioned so far use as starting point the limit of infinitely heavy 
nucleons; corrections due to the finite nucleon mass are then included as a 
power expansion in 1/jVf/v- However, it was already observed in the 70s [40], 
based on calculations using Gaussian wave-functions, that this way one may 
miss important terms. This was further investigated in Ref. [15], where the 
analysis was done model independently and the appearance of these additional 
contributions was related to the Pauli principle in the intermediate AW state, 
while the pion is in flight. This will be discussed in detail in the next section. 
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3 Role of ttNN cuts 



Let us investigate the diagrams of Fig. 3 from the point of view of their irNN 
cut. Then we may write for the corresponding matrix elements 

[ dK„K*dP R P* f(K 2 ,P 2 R ) 

1«NN(Q)-J (2yr)6 Q. K 2 /{2mn) _p2 R/MN+l0 > W 

where the function / contains the reaction specific parts like vertex functions, 
wave functions and transition operators. The only part spelled out explicitly 
is the ttNN propagator, where for simplicity non-relativistic kinematics was 
chosen. Here Q denotes the excess energy with respect to the irNN thresh- 
old. The goal of this study is to compare the full expression given in Eq. 3 
with the corresponding one that emerges when static nucleons are used. This 
corresponds to taking the limit Mn — > oo prior to integration and we may 
write 

^static) {0) = f dK n K^dP R P 2 f{Kl,P 2 R ) _ / m 



" NN {Q) ~ I VvO 6 Q-Klf{2m,)+iQ +0 \M N ) ' (4) 

It was generally assumed that the difference between the integrals of Eqs. (3) 
and (4) can be accounted for in a polynomial expansion in iti^/Mn ■ However, 
this is in general not the case: at pion production threshold I n NN(Q) has 
a branch point singularity. As we will see this leads to a contribution non- 
analytic in iti^/Mn, even below pion threshold. To see this, let us study 
in detail the cut contribution by replacing the ttNN propagator by its delta 
function piece 

(Q - K 2 J{2m n ) - P 2 /M N + iO)- 1 - -iir6(Q - K 2 J{2m^) - P 2 R /M N ) . 
Then we may write 

C%(Q) = -^j ^^f(K n (Q,P 2 n) 2 ,P 2 R )K n (Q,P 2 R) ■ 

for the full contribution, where K™(Q,P%) = ^2m v {Q - P R /M N ) and 
llff atiC \Q) = -i™ n K™(Q) J ^ff(K™(Q) 2 ,P 2 ) 



for the static one, with K° n {Q)\ static = ^2m n Q . At this point we see already 
one important difference between the static and the full treatment: while 
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the imaginary part of the former scales in a completely wrong way, namely 
according to two-body phase-space, the latter shows the proper scaling as 
three-body phase space. 

But this is not all. Also below the 7rNN threshold the static contribution 
gives wrong results. To see this let us focus on the contribution at irNN 
threshold , Q = 0. Although for this kinematics both the above integrals are 
real, the presence of the irNN cut still plays a significant role. To evaluate the 
relevant integral the on-shell momentum K° n (Q,P^) needs to be continued 
analytically to imaginary values using the prescription 



Thus, taking the Mjy — ► oo limit prior to integration is in general not allowed 
in the presence of few-body cuts. 

The natural question is: when does this matter? As explained the men- 
tioned effect originates from the opening of the physical irNN threshold. How- 
ever, this threshold can only matter, if the irNN state is allowed by selection 
rules. In the isospin limit the two nucleons are identical particles that are to 
obey the Pauli principle. It therefore depends on the operator that acts on 
the deuteron wave function to produce the intermediate pion, whether or not 
the NN state is allowed and, consequently, whether the above contributions 
matter. 

Let us first look at ird scattering. The leading operator that contributes 
to the ttN — ► ttN transition is the so-called Weinberg-Tomozawa term oc 
e abc T c , which is a vector in isospin space, but spin and momentum independent. 
Therefore, this operator acting on the deuteron (isospin and spin 1), leads 
to an NN state that is isospin 1 and spin 1, predominantly in an s-wave due 
to the momentum independence of the transition operator. This NN state is 
forbidden by the Pauli principle and therefore all said in the first part of this 
section does not matter and the static approximation gives a good description 
of the leading few-body correction. The same holds, e.g., for ird — > jNN. 

However, there are reactions where we expect the above terms to become 
significant. One example is the reaction — > ir + nn that will be discussed 




With this we get 




whereas 



r (cut,static) 
tvNN 



(0) = . 
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in more detail in the next section. Here the operator acting on the deuteron 
wave function in leading order is the so-called Kroll-Ruderman term, which 
is a vector in both isospin as well as spin space. Thus, a transition to an 
NN pair in the 1 So isovector state, which is allowed by the Pauli principle, 
is possible. This reaction was studied in detail in Ref. [41], and indeed the 
pattern sketched above on general grounds was observed. However, for a 
Pauli allowed intermediate state, the two nucleons will interact. It was found 
that the inclusion of the two-nucleon intermediate state gives a significant 
contribution, however, numerically smaller than the static exchange itself. 
This is different to the case of ird scattering with an isoscalar irN interaction 
that also leads to a Pauli allowed intermediate state. In this case the inclusion 
of the NN interaction at threshold numerically restored the contribution of 
the static exchange [42]. 

4 The reaction — > 7T + nn 

The reaction jd — ► i: + nn was studied intensively already in the 70s — for a 
review see Ref. [43]. At this time only diagrams (al) and (a2) of Fig. 4 were 
included. In Ref. [43] there is only one comment to an unpublished work, where 
pion rescattering (diagrams (cl) and (c2) of Fig. 4) was calculated, however, 
in the static approximation. It is stated that the inclusion of this contribution 
destroys the nice agreement of the calculation based on the one-body terms 
only and therefore it will no longer be considered. Based on the discussion 
of the previous section we now understand, why the static pion exchange 
diagram gave a contribution way too large: in a complete calculation it would 
have been largely canceled by the recoil corrections, since the two-nucleon 
state in diagrams (6), (c) and (d) can go on-shell while the pion is in flight. 

The results of our calculation are shown in Fig. 5. At leading order and 
next-to-leading order only one-body terms contribute (c.f. Fig. 4) with their 
strength fixed by the chiral Lagrangian. The relevant jp — > ir + n vertices are 
momentum independent in both cases and therefore their energy dependence 
is identical. At NNLO there is a counter term for the transition jp — ► ir + n and 
the strength of the one-body operator can be adjusted to data [44]. This gives 
a large fraction of the shift in strength when going from NLO to NNLO. In 
addition the amplitude gets energy dependent [45]. Another source of energy 
dependence comes from the few-body corrections as well as higher partial 
waves that start to contribute at this order. As can be seen from the figure, 
the data is described very nicely in the whole low energy region considered. 

It seems as if the few body corrections, when treated properly, only have 
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Figure 4: Diagrams contributing to ~fd — > ir + nn up to order x 5//2 - As before, 
solid lines denote nucleons, and dashed lines pions. The wavy lines denote 
photons. The hatched areas denote the deuteron wave function and the filled 
circles the NN interaction. 

a minor effect on the event rates for 'yd — > ir + nn, however, this is correct only 
for the total cross section. Especially the neutron momentum distributions 
are sensitive to higher order corrections and those are to be understood to 
very high accuracy, in order to make use of this reaction to extract the nn 
scattering length. 

On the level of neutron momentum distributions diagram (al) leads to 
very specific signals due to the so-called quasi-free production. When all 
particles in the final state go forward, the intermediate proton is very near 
on-shell and the diagram gives a large contribution, which decreases quickly, 
however, as we go away from forward kinematics. In addition, the quasi-free 
production favors large relative momenta of the two neutrons. Near threshold 
this is clear, as — in the center of mass system — the spectator neutron keeps 
on going with half the deuteron momentum, whereas the reaction neutron gets 
decelerated to almost at rest through the production process. 

All diagrams with a final-state interaction, on the other hand, give contri- 
butions peaked at small relative nn momenta and almost insensitive to their 
orientation. This is illustrated in Fig. 6, where the differential rate is shown 
for two different angles as a function of the relative nn momentum Pr. The 
dashed line denotes the distribution where Pr is directed along the beam axis, 
whereas for the solid line it is perpendicular to the beam. In the former case 
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Figure 5: Total cross section of the reaction 7c? — > ir + nn at LO (dashed line), 
NLO (dash-dotted line) and x 5 ^ 2 order (solid line) together with experimental 
data from Ref. [46]. 



the distribution shows a clear two-peak structure — the quasi-free production 
shows up at large Pr and the final state peak shows up at small Pr. In the 
latter case, on the other hand, the quasi free contribution has disappeared and 
only the final state interaction piece remains with basically identical strength. 

The hight and shape of the FSI peak is sensitive to the value of a nn — 
the neutron-neutron scattering length. A systematic study revealed that high 
accuracy data on jd — > ir + nn will allow one to extract the nn scattering 
length with an uncertainty of the order of 0.1 fm which is compatible with that 
estimated for the competing reactions, pd — * nnp [47] and ir~d — >• jnn [13]. 

5 NN -> dir 

As sketched in the introduction, for reactions of the type NN — > NNtt a 
simultaneous expansion in the large initial momentum pthr ~ V m Tr^N, that 
also sets the scale for the typical momenta in the loops, and the pion mass 
is compulsory. Before we go into details in discussing a particular reaction 
channel we would like to briefly illustrate the impact of this. In practice this 
means that momenta and pion masses are to be treated independently in the 
power counting. 

To see how this works let us for example estimate the contributions of the 
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Figure 6: Predicted event rate for the reaction 'yd — > ir + nn at an excitation 
energy of 5 MeV as a function of the nn relative momentum Pr for two 
different orientations of it. The region of small values of Pr are dominated 
by the nn final state interaction (FSI) and that of large values of it by the 
quasi-free production (QF). 



loops shown in Fig. 7. For diagram (a) we then estimate 

P ( P \ 3 (i\ 2 (i\ 2 p 4 l_ 

ft\h) \P 2 J \p) (47r) 2 ~ M$j ' 

where the different terms refer to the ttN — * ttN vertex, the three ttNN ver- 
tices, the two pion propagators, the two nucleon propagators, and the integral 
measure, in order. Each individual piece was expressed by the dimension- 
ful parts, where momenta were identified with their typical values. For more 
details we refer to Appendix E of Ref. [24]. To come to the order estimate 
we used 47r/ 7r ~ Mjy and dropped an overall factor of l/f% common to all 
production amplitudes. On the other hand we find for diagram (6) 

j fm n \ 2 1 1 m% \ 1 / p \ ml 

ml (4vr)2 j p 2 \ fj ~ pM% ' 

The expression in the curly bracket refers to the pion loop — it contains two 
irN — * irN vertices, one nucleon and one pion propagator as well as the integral 
measure — however, in this case the typical momentum is of the order of m n 
instead of pthr as in the previous example. The reason is simply that one may 
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Figure 7: Two typical pion loops that contribute to NN — > NNir + . Diagram 
(a) starts to contribute at NLO whereas diagram (6) starts to contribute at 
N 4 LO. 



a) 



b) 



c) 



Figure 8: Tree level diagrams that contribute to pp — ► dir + up to NLO. Solid 
lines denote nucleons, dashed ones pions and the double line the propagation 
of a Delta-isobar. 



choose in the loop the momenta such that the large momentum does not run 
through the pion. Then the large scale does not appear in the loop at all since 
the leading heavy baryon propagator feels energies only. Outside the curly 
bracket the momentum transfer is large and therefore the pion propagator as 
well as the pion vertex appear with p ~ pthr- 

If we compare the two order asignments we observe that under the assump- 
tion that momenta are of order m n , both expressions appear at the same order 
and are therefore expected to be of similar order of magnitude. However, if 
we assume p to be of order p t hr, then diagram (a) is of order m^/M^ and dia- 
gram (6) is of order (m^/Miy) 5 ' 2 , which corresponds to a relative suppression 
of (m^/MTv) 3 / 2 ~ 1/20. An explicit calculation [19] revealed an even larger 
suppression of diagram (6), which turned out to be suppressed by a factor of 
50 compared to (a). For more examples we refer to Ref. [24]. 
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Figure 9: Irreducible pion loops with nucleons only that start to contribute to 
AW -> AWvr at NLO that were considered in Ref. [25]. 

As already discussed in the context of nd scattering, in some cases only 
the sum of various diagrams can give meaningful results, since the individual 
diagrams depend on the choice made for the pion field. Due to the reordering 
of loop diagrams, the members of those invariant subgroups are not neces- 
sarily of the same chiral order anymore, in contrast to the original Weinberg 
counting. That, regardless this, also the new scheme gives meaningful results 
is demontrated in Ref. [48]. 

The tree level amplitudes that contribute to pp — ► di: + are shown in Fig. 8. 
In Ref. [25] all NLO contributions of loops that start to contribute to AW — > 
NNir at NLO were 1 calculated in threshold kinematics — that is neglecting 
the distortions from the AW final- and initial state interaction and putting all 
final states at rest. At threshold only two amplitudes contribute, namely the 
one with the nucleon pair in the final and initial state in isospin 1 (measured, 
e.g., in pp — > ppir°) and the one where the total AW isospin is changed from 1 
to (measured, e.g., in pp — > dir + ). It was found that the sum all loops that 
contain A-excitations vanish in both channels. This was understood, since 
the loops were divergent and at NLO no counter term is allowed by chiral 
symmetry. On the other hand the nucleonic loops were individually finite. It 
was found that the sum of all nucleonic loops that contribute to pp — > ppn 
vanish, whereas the sum of those that contribute to pp — > dir + gave a finite 
answer. The resulting amplitude grows linear with the initial momentum. In 
Ref. [49] it was pointed out that this growth of the amplitude is problematic: 
when evaluated for finite outgoing AW momenta, the transition amplitudes 
turned out to scale as the momentum transfer. Especially, the amplitudes 
then grew linearly with the external AW momenta. As a consequence, once 
convoluted with the AW wave functions, a large sensitivity to those was found, 
in conflict with general requirements from field theory. The solution to this 
puzzle was presented in Ref. [50] and will be reported now. 

'in a scheme with two expansion parameters — here m n and pthr — loops no longer 
contribute at a single order but at all orders higher than where they start to contribute. 
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Figure 10: The ttN — > 7riV transition vertex: definition of kinematic variables 
as used in the text. 

The observation central to the analysis is that the leading ttN — > 7riV 
transition vertex, as it appears in Fig. 8a, is energy dependent. Using the 
notation of Fig. 10 its momentum and energy dependent part may be written 
as [2] 



For simplicity we skipped the isospin part of the amplitude. The first term 
in the last line denotes the transition in on-shell kinematics, the second the 
inverse of the outgoing nucleon propagator and third the inverse of the incom- 
ing nucleon propagator. First of all we observe that for on-shell incoming and 
outgoing nucleons, the the ttN — > irN transition vertex takes its on-shell value 
2m T — even if the incoming pion is off-shell, as it is for diagram (a) of Fig. 
8. This is in contrast to standard phenomenological treatments [51], where lo 
was identified with 777.^/2 — the energy transfer in on-shell kinematics — and 
the recoil terms were not considered. Note, since p^ hr /M^ = m w the recoil 
terms are to be kept. 

A second consequence of Eq. (5) is even more interesting: when the 
ttN — > ttN vertex gets convoluted with NN wave functions, only the first 
term leads to a reducible diagram. The second and third term, however, lead 
to irreducible contributions, since one of the nucleon propagators gets can- 
celed. This is illustrated in Fig. 11, where those induced topologies are shown 
that appear, when one of the nucleon propagators is canceled (marked by the 
filled box) in the convolution of typical diagrams of the NN potential with the 
NN — > NNtt transition operator. Power counting gives that diagram (6) and 
(c) appear only at order N 4 LO and N 3 LO, respectively. However, diagram (a) 
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a) b) c) 

Figure 11: Induced irreducible topologies, when the off-shell terms of Eq. (5) 
hit the NN potential in the final state. The filled box on the nucleon line 
denotes the propagator canceled by the off-shell part of the vertex. 

starts to contribute at NLO and it was found in Ref. [50] that those induced 
irreducible contributions cancel the finite remainder of the NLO loops in the 
pp — > dir + channel. Thus, up to NLO only the diagrams of Fig. 8 contribute 
to pp — > dir + , with the rule that the irN — > irN vertex is put on-shell. 

The result found in Ref. [50] is shown in Fig. 12 as the solid line, where 
the total cross section (normalized by the energy dependence of phase space) 
is plotted against the normalized pion momentum. The dashed line is the 
result of the model by Koltun and Reitan [51], as described above. The data 
sets are from TRIUMF [52], IUCF [53], and COSY [54]. 

6 Comparison to phenomenological works 

All recent phenomenological calculations for NN — > NNir add additional di- 
agrams to the model of Ref. [51]. Here we will focus only on pp — ► dn + . 
Phenomenological calculations for this reaction in near threshold kinematics 
are given, e.g., in Ref. [55] and Ref. [56]. In both works in addition to the 
diagrams of Ref. [51] some A-loops as well as additional short range contri- 
butions are included — heavy meson exchanges for the former and off-shell 
ttN scattering 2 for the latter. Based on this the cross section for pp — > dir + is 
overestimated near threshold. How can we interpret this discrepancy in light 
of the discussion above? 

First of all, the NLO parts of the A-loops cancel, as was shown already 
in Ref. [25]. However, in both Refs. [55,56] only one of these diagrams was 
included and, especially for Ref. [55], gave a significant contribution. In 
addition, in the effective field theory short ranged operators start to contribute 
only at N 2 LO. The only diagram of those NLO loops shown in Fig. 9 that is 

2 That those are also short range contributions is discussed in Ref. [24]. 
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Figure 12: Comparison of the results of Ref. [50] to experimental data for 
NN — > dir. The dashed line corresponds to the model of Koltun and Rei- 
tan [51], whereas the solid line is the result of the ChPT calculation of Ref. [50]. 
The estimated theoretical uncertainty (see text) is illustrated by the narrow 
box. The data is from Refs. [52] (open circles), [53] (filled circles) and [54] 
(filled squares). The first data set shows twice the cross section for pn — > dn 
and the other two the cross section for pp — > dir + . 

effectively included in Ref. [56] is the fourth, since the pion loop there can 
be regarded as part of the irN — > irN transition T-matrix. However, as 
described, the contribution of this diagram gets canceled by the others shown 
in Fig. 9 and the induced irreducible pieces described above. Therefore, the 
physics that enhances the cross section compared to the work of Ref. [51] in 
Refs. [55,56] is completely different to that of Ref. [50]. However, only the 
last one is field theoretically consistent as explained in the previous section. 

The natural question that arises is that for observable consequences. As 
explained, in the effective field theory calculation the near threshold cross sec- 
tion for pp — > dir + is basically given by a long-ranged pion exchange diagram, 
whereas the phenomenological calculations rely on short ranged operators with 
respect to the AW system. Obviously those observables are sensitive to this 
difference that get prominent contributions from higher partial waves in the 
final iViV system. We therefore need to look at the reaction pp — > pnir + . Un- 
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fortunately, the total cross section for this reaction is largely saturated by NN 
5-waves in the final state (see, e.g., Fig. 17 in Ref. [24]). On the other hand, 
linear combinations of double polarization observables allow one to remove the 
prominent components and the sub-leading amplitudes should be visible. We 
therefore expect from the above considerations that the phenomenological cal- 
culations give good results for polarization observables for pp — > dir + , whereas 
there should be deviations for some of those for pp — > pmr + . Predictions for 
these observables were presented in Ref. [57] and indeed the ir + observables 
with the deuteron in the final state, reported in Ref. [58], are described well 
whereas there are discrepancies for the pn final state (see Fig. 24 of Ref. [24] ) . 
For the corresponding data see Ref. [59]. 

It remains to be seen how well the same data can be described in the 
effective field theory framework. Up to NNLO the number of counter terms is 
quite low: there are two counter terms for pion s-waves, that can be arranged 
to contribute to pp — > ppn° and pp — > d,7r + individually, and then there is 
one counter term for pion p-waves, that contributes only to a small amplitude 
in charged pion production [23]. On the other hand there is a huge amount 
of even double polarized data available [58-60] — and there is more to come 
especially for pn — > ppir~ [61]. 

7 Dispersive corrections to a n d 

Let us come back to ird scattering at threshold. The corresponding scattering 
length was presented in Eq. (2). In the introductory sections we exclusively 
focused on the real part. However, now we are in the position to also discuss 
the imaginary part, which is closely linked to the reaction NN — ► dir through 
unitarity and detailed balance. One may write 

47rIm(o ff d) = lim q {a(iTd -► NN) + a(ird -► ~fNN)} , (6) 
<2^0 

where q denotes the relative momentum of the initial ird pair. The ratio 
R = lim 9 _ ((r(ird -> NN)/a(Trd -> jNN)) was measured to be 2.83 ± 0.04 
[62]. At low energies diagrams that lead to a sizable imaginary part of some 
amplitude are expected to also contribute significantly to its real part. Those 
contributions are called dispersive corrections. As a first estimate Bruckner 
speculated that the real and imaginary part of these contributions should be of 
the same order of magnitude [63]. This expectation was confirmed within Fad- 
deev calculations in Refs. [64]. Given the high accuracy of the measurement 
and the size of the imaginary part of the scattering length, another critical 
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Figure 13: Dispersive corrections to the nd scattering length. 

look at this result is called for as already stressed in Refs. [65,66]. A con- 
sistent calculation is only possible within a well defined effective field theory 
- the first calculation of this kind was presented in Ref. [67] and is briefly 
sketched here. 

To identify the diagrams that are to contribute we first need to specify 
what we mean by a dispersive correction. We define dispersive corrections 
as contributions from diagrams with an intermediate state that contains only 
nucleons, photons and at most real pions. Therefore, all the diagrams shown 
in Fig. 13 are included in our work. On the other hand, all diagrams that, e.g., 
have Delta excitations in the intermediate state do not qualify as dispersive 
corrections, although they might give significant contributions [31]. 

The hatched blocks in the diagrams of Fig. 13 refer to the relevant tran- 
sition operators for the reaction AW — * NNir depicted in Fig. 8. Also in 
the kinematics of relevance here the irN —* irN transitions are to be taken 
with their on-shell value 2m n . Using the CD-Bonn potential [68] for the AW 
distortions we found for the dispersive correction from the purely hadronic 
transition 

a^Y = (-6.3 + 2 + 3.1 - 0.4) x 1(T 3 m' 1 = -1.6 x 1(T 3 m" 1 , (7) 

where the numbers in the first bracket are the individual results for the dia- 
grams shown in Fig. 13, in order. Note that the diagrams with intermediate 
AW interactions and the crossed ones (diagram (c) and (d)), neither of them 
included in most of the previous calculations, give significant contributions. 
The latter finding might come as a surprise on the first glance, however, please 
recall that in the chiral limit all four diagrams of Fig. 13 are kinematically iden- 
tical and chiral perturbation theory is a systematic expansion around exactly 
this point. Thus, as a result we find that the dispersive corrections to the 
7tg? scattering length are of the order of 6 % of the real part of the scatter- 
ing length. This number is fully in line with the expectations from power 
counting, which predicted a relative suppression of the dispersive corrections 
compared to the leading double scattering term — diagram (6) of Fig. 1 - 
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of the order of (m n /MN) 3 / 2 ~ 5 %. Note that the same calculation gave very 
nice agreement for the corresponding imaginary part [67]. 

In Ref. [67] also the electro-magnetic contribution to the dispersive cor- 
rection was calculated. It turned out that the contribution to the real part 

was tiny 0.1 x lO^m^ 1 — while the sizable experimental value for the 

imaginary part (c.f. Eqs. (2) and 6) was described well. 

To get a reliable estimate of the uncertainty of the calculation just pre- 
sented a NNLO calculation is necessary. At that order a counter term appears 
for pions at rest that can be fixed from NN — > NNtt, as indicated above. 
For now we can only present a conservative estimate for the uncertainty by 
using the uncertainty of order 2m 7r /M/v one has for, e.g., the sum of all direct 
diagrams to derive a Aa^ p of around 1.4 x lO^m^ 1 , which corresponds to 

about 6% of Re (a^ 3 ^ . However, given that the operators that contribute 
to both direct and crossed diagrams are almost the same and that part of the 
mentioned cancellations is a direct consequence of kinematics, this number for 
Aa^ p is probably too large. 

In Ref. [67] a detailed comparison to previous works is given. Differences 
in the values found for the dispersive corrections were traced to the incomplete 
sets of diagrams included in those phenomenological studies. 

8 Summary and Outlook for pion reactions 

In the lectures recent progress in our understanding of elastic and inelastic 
pion reactions on the two-nucleon system was presented. The central reaction 
discussed was nd scattering at threshold. Arguments were given that in the 
years to come one should be able to calculate the nd scattering length with 
sufficient accuracy to use the reaction as one of the prime sources for the 
isoscalar scattering length a + . To reach this not only significant progress was 
necessary for the coherent ltd scattering but also the reactions NN — > NNtt 
need to understood. In the future also the role of isospin violation on ird 
scattering needs to be investigated further as stressed in Ref. [35]. 

The process NN — > NNtt is a puzzle already since more than a decade. 
Given the progress presented above we have now reason to believe that this 
puzzle will be solved soon. This mentioned results could only be found, because 
a consistent effective field theory was used. For example, the potential problem 
with the transition operators of Ref. [25], pointed at in Ref. [49], would always 
be hidden in phenomenological calculations, since the form factors routinely 
used there always lead to finite, well behaved amplitudes. The very large 
number of observables available for the reactions NN — ► NNtt will provide a 
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Figure 14: Comparison of different variants for the AiV interaction to the 
available data at low energies. The solid line, dashed line, and the light area 
are the results of Refs. [77-79], in order. The data are from Refs. [75, 76]. 

non-trivial test to the approach described. 

Once the scheme is established, the same field theory can be used to analyze 
the isospin violating observables measured in pn — > dir° [69] and dd — > air [70] . 
First steps in this direction were already done in Ref. [71] for the former and 
in Refs. [72, 73] for the latter. 

In the lectures also the reaction jd — ► Tr + nn was discussed. Not only gave 
those studies a further confirmation that we understand the few body dynam- 
ics well within ChPT, but it also promises to become an ideal reaction for the 
extraction of the nn scattering length with high accuracy. The corresponding 
measurements could be performed at HIGS [74]. 

9 Some remarks on strangeness production 

As described in the previous sections a lot is already known about the proper- 
ties and dynamics of systems composed of light quarks only. However, much 
less is known about the scattering of systems with strangeness, especially for 
low energies. The reason is of experimental nature: the lifetime of particles 
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with strangeness is typically too short to allow for secondary scatterings at 
low energies, necessary to get information on low energy scattering. 

Here we will focus on the hyperon-nucleon (YN) system. The poor status 
of our information on the YN interaction is most obviously reflected in the 
present knowledge of the A7V scattering lengths. Attempts in the 1960's to 
pin down the low energy parameters for the 5-waves led to results that were 
afflicted by rather large uncertainties [75, 76]. In Ref. [76] the following values 
are given for the singlet scattering length a s and the triplet scattering length 
a t 

a* = -1-8 fm and a* = -1.6 j^'g fm, (8) 

where the errors are strongly correlated. The situation of the corresponding 
effective ranges is even worse: for both spin states values between and 16 
fm are allowed by the data. Later, the application of microscopic models for 
the extrapolation of the data to the threshold, was hardly more successful to 
pin down the low energy parameters. For example, in Ref. [77] one can find 
six different models that equally well describe the available data but whose 
(iS-wave) scattering lengths range from -0.7 to -2.6 fm in the singlet channel 
and from -1.7 to -2.15 fm in the triplet channel. To illustrate this point in Fig. 
14 we show a comparison of model f of Ref. [77] (dark solid line), the Jiilich 
'04 model [78] (dashed curve), and the result from the recent effective field 
theory approach of Ref. [79] to the world data. 

The natural alternative to scattering experiments are production reactions. 
However, the central insights of the previous sections were that only within a 
consistent field theory reliable calculations can be performed for the reactions 
under considerations. On the other hand, as stressed in the introduction, any 
strangeness production reaction involves momenta that do not allow for an 
expansion along the lines just discussed, since the corresponding expansion 
parameter in this case would be larger than 1/2. Does this mean that one can 
learn nothing from a study of strangeness production off two nucleon systems? 

Not at all. For one thing, an investigation of baryon and meson resonances 
does not need any detailed knowledge on the production mechanism — most 
of the relevant information is contained in the Dalitz plots — and the com- 
ment of the previous paragraph applies only to this part. But one can learn 
even more from the production reactions by using that the momentum transfer 
in those reactions are large. Since this leads to an effectively point-like pro- 
duction operator, one may employ dispersion theory to relate invariant mass 
spectra of production reactions to elastic scattering data. Obviously, for this 
no knowledge on the production operator is necessary whatsoever. Then, in 
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contrast to above, one works with the typical outgoing relative momentum in 
units of the momentum transfer as expansion parameter. 

In the remainder of this text I will only focus on how to extract scatter- 
ing parameters from production reactions. One aspect of spectroscopy, also 
discussed in the lectures, namely that of scalar mesons, was described in the re- 
cent conference proceeding [80] in very much detail and it will not be repeated 
here. 

The use of dispersion theory was very common in the 50s and the basis for 
the study to be described now was worked out already then [81]. A controlled 
method of extraction of scattering lengths from production reactions opens up 
the opportunity to measure scattering parameters also for unstable states. As 
an example we will discuss the option to measure the AN scattering lengths 
from NN -> KAN and -> KAN. In Refs. [82-84] it was shown, how to 
derive an integral representation of the scattering length of a pair of outgoing 
particles (here we show the formula relevant for neutral final states, as is 
relevant for AN; in the presence of Coulomb interactions the equation has to 
be modified — see Ref. [83]): 




where as denotes the spin cross section for the production of a A-nucleon 
pair with invariant mass w! 2 — corresponding to a relative momentum p' — and 
total spin S. In addition t = (p\ — Pk+) 2 i w hh p\ being the beam momentum, 
777.Q = (ttt-a + mjv) 2 , where m\ (tun) denotes the mass of the Lambda hyperon 
(nucleon), and m max is some suitably chosen cutoff in the mass integration. 
In Ref. [82] it was shown that it is sufficient to include relative energies of the 
final AiV system of at most 40 MeV in the range of integration to get accurate 
results. P denotes that the principal value of the integral is to be used and 
the limit has to be taken from above. 

The formula as given is applicable if there are no significant effects from 
crossed channels — this can be monitored by a Dalitz plot analysis — and 
only a single partial wave contributes. With respect to the angular momen- 
tum this can be achieved by a proper choice of m max . In order to select a 
single spin state for the outgoing two-particle system polarization observables 
are necessary. In Refs. [82, 84] the relevant observables are identified for the 
reactions NN — > NKY and jd —> KAN, respectively. The former class of 
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Figure 15: Illustration of the uncertainties of extraction of scattering parame- 
ters from scattering data (left panel) , where an extrapolation is necessary, and 
from production data (right panel), where the data needs to be interpolated. 
The data was taken from Refs. [76, 85] for the right panel and from Ref. [86] 
for the left one. 

reactions can be measured at COSY [61, 87] and the latter either at J-Lab [88], 
MAMI [89], or ELS A [90]. 

The essential advantages of the use of Eq. (9) to extract the scattering 
lengths of unstable particles compared to a determination from scattering 
experiments are, e.g., for AiV scattering: 

• Instead of an extrapolation of data, the scattering length is found from 
an interpolation of an invariant mass spectrum, which is theoretically 
much better controlled. This is illustrated in Fig. 15. 

• The integral representation gives a result for the scattering length with- 
out any assumption on the energy dependence of the hyperon-nucleon 
interaction. This opens the possibility to fix the scattering lengths from 
production reactions and then use scattering data to fix, e.g., the effec- 
tive range. 

• Since Eq. (9) is derived from a dispersion integral, a controlled error 
estimate is possible. A systematic study revealed that for the kinemat- 
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ics relevant for the production of the AN system at low energies, the 
uncertainty of Eq. (9) was found to be 0.5 fm [82,84]. Sources of this 
uncertainty are the possible energy dependence of the production op- 
erator, the influence of the upper limit of integration, and the possible 
influence of crossed channel effects. 

10 Summary 

In these proceedings various reactions on few nucleon systems were discussed. 
It was demonstrated that due to significant advances in the technologies of 
effective field theories, high precision calculations became possible for hadronic 
reactions even on few-nucleon systems. At the same time a clear connection 
to QCD is provided. 

For the production of more heavy systems, like those that contain the 
strangeness degree of freedom, off two nucleons no effective field theory is 
developed yet and therefore many reactions are still being analyzed using 
models. However, also for this class of reactions some aspects can be analyzed 
in a model independent way. The AN scattering lengths were discussed as an 
example. 

We are now in a phase, were both theory and experiment are advanced 
sufficiently that we should understand much better how QCD influences low 
energy nuclear dynamics in the upcoming years. Especially the symmetry 
breaking sector — violation of isospin as well as violation of the flavor SU(3) 
- promises deep insights into the mechanisms of strong interactions. 
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